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Abstract 

The Ml decay widths of charm and bottom mesons and their excited 
states are sensitive to the relativistic aspect of the quark current operators 
and obtain a significant contribution from the interaction current that is 
associated with the scalar confining interaction between heavy quarks (Q) 
and light antiquarks (q). Consequently they provide direct evidence on 
the latter. The spectra and the Ml transition widths of the D, Ds and 
B, Bs mesons and their orbital excitations are calculated here within the 
framework of the Blankenbecler-Sugar equation, which allows a covari- 
ant treatment, while retaining a formal similarity to the nonrelativistic 
approach. The hyperfine interaction in the Qq mesons is modeled by 
screened gluon exchange, which shares many features with the instanton 
induced interaction. The quality of the model is tested by calculation of 
the spectra and Ml decay widths of charmonium. 



* talahdeOrock. helsinki . fi 
^ nyfalt @rock. helsinki . fi 
■f riska@science.helsinki.fi 



1 



1 Introduction 



A dynamical description of the spectrum of the heavy flavor mesons Qq has 
a priori to be relativistic as the confined light quarks and antiquarks have ve- 
locities close to that of light. Combination of this requirement, with the non- 
relativistic instantaneous aspect of the effective linear [Q (or near linear 
confining interaction suggests that a 3-dimensional quasipotential reduction of 
the Bethe-Salpeter equation should provide the most obvious framework for a 
realistic description. This is a fortiori the case as the Bethe-Salpeter equation 
with instantaneous interactions, while commonly employed, fails to provide 
satisfactory results for the Ml transitions of the analogous heavy quarkonia Q. 

Among the many choices of quasipotential equations the Blankenbecler- 
Sugar equation H, |^ suggests itself because of its transparent formal similarity 
to the nonrelativistic Lippmann-Schwinger equation, to which it reduces in the 
adiabatic limit. Moreover it admits direct employment of the conventional (rel- 
ativistic) models for the hyperfine interaction between quarks and antiquarks, 
as well as the linear confining interaction along with its obvious relativistic 
modifications. The key advantage of the Blankenbecler-Sugar equation over the 
nonrelativistic or "relativized" versions of the nonrelativistic wave equation is 
the unapproximated treatment of the relativistic two-particle kinetic energy op- 
erator of the Qq system. The relativistic modifications to the interaction are in 
comparison modest, and may be treated in an approximate fashion, when conve- 
nient. This feature is shared with the Gross reduction of the Bethe-Salpeter 
equation, which has been applied to heavy flavor mesons in ref. Q . 

Given this situation, several different versions of the Blankenbecler-Sugar 
equation have been employed in the literature We shall here use the sym- 
metrical version of ref. to calculate both the spectra and the Ml decay 
widths of the D and Dg as well as the B and Bg mesons and their orbital exci- 
tations. The formalism then differs somewhat from that used to describe the Qq 
system in refs. . The main advantage of the version used here is that it leads 
to a Lippmann-Schwinger type wave-equation. It differs most obviously from 
that used in r efs. [ill in the appearance of the "minimal relativity" square 
root factors y^m/E in the effective interaction potential in the wave equation. 

The radiative transitions of the ground state Qq mesons may be described 
phenomenologically by effective field theory methods jl^, ^ A complete 
dynamical model is however required for any comprehensive description of the 
decays of the full spectrum of Qq states. The static quark model is insuffi- 
cient here, as it is expected to overpredict the width of the Ml transition for 
Z)+* by a large factor in analogy to its large overprediction of the width 

for the Ml decay J/ip — > 7yc7, and which is resolved only by employment of the 
relativistic magnetic moment operator for the heavy quarks along with the inter- 
action current operator that is associated with the scalar confining interaction 

In comparison to heavy quarkonia the heavy- light Qq mesons are more com- 
plex systems in view of the still unsettled form of the hyperfine interaction for 
the light quarks and concomitantly their empirically still only poorly known 
spectrum. While the branching ratios of the Ml decays of the charmed vector 
mesons are known, their still unknown total decay widths prevent an absolute 
empirical determination of their Ml decay widths. In spite of - and because of 
- this situation it remains a challenge to decode the structure and the dynamics 
of these systems, as this is a main key to the understanding of the nature of the 



2 



confining interaction and the hyperfine interaction between constituent quarks. 

There are good theoretical |jl|, |l^ and phenomenological reasons to ex- 
pect the confining interaction to be Hnear with a scalar spinor structure, al- 
though mixtures of combinations of scalar and vector coupling terms continue 
to be employed |l^ . For the heavy light Qq systems the hyperfine interaction 
is likely to be formed as a combination of the screened perturbative gluon ex- 
change interaction and the direct instanton induced interaction which 
share many qualitative features. The latter appears in the instanton liquid 
model description of the QCD vacuum which is supported by numerical 
lattice calculations |22) . In the case of heavy light quarkonia the relativis- 
tic corrections to the gluon exchange interaction cause an appreciable damping 
of that interaction, besides the screening of the the color hyperfine coupling 
strength. 

We here consider the Ml decay rates of the Qq systems, because these are 
expected to be overpredicted by the static quark model, and are likely to pro- 
vide decisive information on the nature of the confining interaction. Because 
the constituent quark masses appear in different combinations in the magnetic 
transition operators for the different charge states of the Qq mesons the Ml tran- 
sition rates should also provide information on the constituent quark masses. 
We show that the decay rates predicted by the static quark model for the Ml 
transition D* Dj is strongly reduced with (1) employment of the relativistic 
magnetic moment operator for the quarks in combination with (2) the interac- 
tion current, that appears if the confining interaction couples as a scalar. This 
result falls in line with that for the corresponding Ml decay rate problem of 
heavy quarkonia [|l6|. For qq systems (i.e. systems with equal quark and anti- 
quark masses) neither the gluon exchange interaction nor the instanton induced 
interaction gives rise two any two-body current operator with a spin-fiip part. 
The spin-flip component of those exchange current operators that appear in the 
case of Qq systems usually have much smaller matrix elements compared to 
the matrix elements of the two-body current that is associated with the scalar 
confining interaction. 

The effective interaction potential between the heavy quark and the light 
antiquark is accordingly described by a scalar linear confining -I- one-gluon ex- 
change interaction with a hyperfine interaction formed of the perturbative gluon 
exchange interaction. The parameters of this model are, to the extent possible, 
fixed to values obtained by lattice calculations as well as to the known part of the 
D and Dg meson spectra, and then appropriately scaled to the bottom mesons. 
Finally the quality of the model is tested by a calculation of the spectrum and 
the Ml decay widths of charmonium. 

With the framework of the Blankenbecler-Sugar equation and the Qq inter- 
action models above it is possible to obtain spectra for both charmonium and 
the heavy-light mesons, which agree overall with the empirical spectra, where 
these are known. Similarly the description of the (known) Ml decay rates of 
the J/ip and the ip' are well described. The Ml decay rates for the D and B 
mesons are predictions, as absolute empirical determination of these are yet to 
be made. The present result that it is possible to achieve a satisfactory descrip- 
tion of both the spectra and Ml decay rates of heavy quarkonia and heavy light 
mesons differs from that of ref. . The key point here is the employment of an 
unapproximated quark current operator. 

This report falls into 7 sections. Section ^ contains a description of the 
application of the Blankenbecler-Sugar equation to the Qq and QQ systems. 
Section |^ contains a description of the interaction parts of the Hamiltonian and 
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the current operators. The calculated D,Ds and B,Bs meson as well as the 
cc spectra are described in section ^. Section ^ describes the quark current 
operators, including the two-quark current that arises along with a linear scalar 
confining interaction. The calculated Ml transition rates are reported in section 
^. Section ^ contains a concluding discussion. 



The Blankenbecler— Sugar equation applied to 

Qq systems 



The Blankenbecler-Sugar equation for a Qq system may be expressed as an 
eigenvalue equation of the form 

+v]il: = e^. (1) 
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Here rrir is the reduced mass of the heavy quark [Q) mass M and the light 
the antiquark [q) mass m: (m^ — Mm/{M + to)). The the eigenvalue e is 
related to the energy E of the Qq meson system as 

^ [E^ -{M-mf][E^ -[M + mf] 
^ SnirE^ ■ ^ 

The interaction operator V, which need not be a local operator, may be obtained 
from the Qq irreducible quasipotential V (in momentum space) as p3; 



Af + TO 



Here the function W is defined as 

W{p)^EQip) + E^ip), (4) 



with Eq{p) — \J + p ^ and Eq(j)) ~ \J m? + p ^ respectively. In the Born 
approximation the quasipotential V equals the Qq invariant scattering amplitude 
T, and thus a constructive relation to field theory obtains. 

Numerical solution of the wave equation (|l|) yields a value for the eigenvalue 
e. The energy (rest mass) of the Qq meson is then obtained as 



E ^ JnP + m'^ + Awirt + 2^JNPm? + 2mr{A'P + 'm^)e + Amle^. (5) 



In the equal mass case M = m this expression reduces to i? = ■\/4to(to + e). 

The relation between the (total) energy E of the Qq system and the eigen- 
value e in the equation (1) shows that the relativistic treatment of the kinetic 
energy term in the Hamiltonian leads to a lowering of the former in comparison 
to what a nonrelativistic treatment would have implied. This is a consequence 
of the effective weakening of the repulsive kinetic energy operator, that results 
when the quadratic terms are replaced by square roots of p^ . Despite the for- 
mal similarity between the Blankenbecler-Sugar equation and the Schrodinger 
equation, the former employs a quadratic mass operator. 
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The main component of the interaction between the Q and q quarks is the 
hnear confining interaction, which is well defined only in the position space 
representation, which we shall therefore adhere to. While the Blankenbecler- 
Sugar equation takes full account of the relativistic two-quark kinetic energy 
operator, the relativistic aspects of the hyperfine interaction operators will be 
treated perturbatively. This treatment is more than adequate for the heavy 
quark, but less so for the light one. It does, however, allow retention of the 
conventional operator structure of the interaction operator between the quarks. 



3 The Interaction Hamiltonian for Qq systems 



3.1 The confining interaction 

The linear confining interaction will here - on the basis of compelling indications 
- be assumed to be a scalar in the spinor representation |ll6| , ^| . To second order 
in the inverse quark masses this interaction, as calculated from the invariant 
amplitude according to (||) , takes the form 

^ '^^\ 2 ml) ^Mmr r AW^nfi 

r SM^'m,'' 

Here c is the string tension parameter, the value of which is ^ 1 GeV/fm. 
The mass coefhcient m2 is defined as 



1 _ VM^ + + Mm 
"12 ^/3Mm 



(7) 



The spin-independent term that is proportional to 1/r is a consequence of the 
square root factors in (^). Without those factors the factor 3/2 in the momentum 
dependent term in the first bracket on the r.h.s, of would be 1. The terms of 
second order in the inverse quark masses in (^) arc implied by scalar coupling 
for the confining interaction. As the antisymmetric spin-orbit interaction has 
no diagonal matrix elements between any of the states in the S- and P-shells, 
and only connects spin singlet and triplet states with J = L it contributes to 
the energy only at quartic order in the inverse quark masses. 

The term of second order in p in the spin-independent term in (^) is un- 
realistically large for light quarks, but is substantially moderated by inclusion 
of the next term in the (asymptotic) (w/c)^ expansion. Similarly the spin-orbit 
interaction in is unrealistically large for light quarks. Both of these explicitly 
velocity dependent terms are strongly moderated by the higher order terms in 
the (asymptotic) expansion in u^. This is seen directly by comparison to the 
corresponding terms in the the quartic correction term to (^, which takes the 
form 



16mf^ r 16to^jj r 32m|^ 
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where Q12 is the quadratic spin-orbit interaction operator 
Q12 = {Sq ■ Ldq ■ L + Oq ■ Luq ■ L)/2. 



(8) 



(9) 



The quartic component of the antisymmetric spin-orbit interaction has been 
dropped here, as it contributes to the energies of the meson states only to sixth 
order. Above the mass coefficients rriAj are defined as 



1 


'{M + nif 






1 


'?,{M + mf 






1 


■5(M + m)4 






1 


■ (M + m)4 







8{M + mf 



16(M + m)2 



3(M + 



1 



(10) 

(11) 
(12) 
(13) 



The scalar structure of the confining interaction implies functional relations 
between the different spin components of the potential. Thus e.g. in the static 
limit the spin-orbit component may be calculated from the central component 
as 

VLs{r) = --(_ + _) -—v,{r). (14) 



4 V 
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Here Vc{r) is the central confining potential and WLs(r) is the coefficient function 
for the spin-orbit coupling operator S ■ L. When the static limit is not invoked, 
the corresponding relation is more complicated. To a good approximation it 
may however be cast in the form 



VLs{r) 



2c d_ 

nr dr 



dr'r'^ 



dkk'^jo{kr)jo{kr') 



{eg + MY (eg + m)2 



(15) 



where eg = -f k'^ /A and Cq = m? -{- k^/A. This relation implies that the 

static approximation to the spin-orbit interaction in (^) represents a considerable 
overestimate of the net spin-orbit interaction in the case of light constituent 
quarks. 

A similar more accurate representation of the quadratic spin-orbit compo- 
nent of the confining interaction obtains with the replacement 



1 



7Q1 



2c_9 i_a -"^ 

nr dr r dr 



dr'r'^ 







dkk 



1 jo{kr)jo{kr') 
(eQ + M)2(e,- + m)2 



(16) 
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The expansion of the interaction is an asymptotic series, which has to 
be truncated. As in the case of constituent quarks the successive terms have 
increasing amphtudes of alternating signs, this series is but a poor representation 
of the full spinorial structure of the interaction, the velocity dependence of which 
is well behaved. We shall therefore truncate the interaction to second order 
and retain only the local spin-orbit and quadratic spin-orbit interaction terms 
of quartic order. The latter is numerically insignificant for the low angular 
momentum states considered here. 

To this order the most realistic treatment of the velocity dependent correc- 
tion —3/2^2/777,2 to the linear potential cr in ^ would be to take it into account 
as a mass shift, 

m* ^ m* + -cr{ — f, (17) 

2 TO2 

in the kinetic term in (^. We shall treat this reduction of the kinetic term 
phenonienologically by subtraction of a constant (6) from the kinetic term. The 
value of this constant is determined by a fit of the calculated spectra to the 
experimental one. Without such a constant the kinetic term in the Hamiltonian 
becomes unrealistically large as becomes evident below. 



3.2 The one-gluon exchange interaction 



In the case of heavy quarkonia the perturbative gluon exchange interaction forms 
an important component of the hyperfine interaction between quarks JTsf . To 
second order in the inverse quark masses (with exception for the quadratic spin- 
orbit interaction, which is quartic) that interaction, after multiplication of the 
square root factors in (||), takes the form: 

4 f 1 37r ^12). . p2 



3 ^ 1 r 2to2 2Mmr 



9 Mm iMmr"^ 
2a, i M"^ +m^ 1 ' 



S-L 



3r3 [ 2M2to2 TV/to 
a, M2 — m^ ,^ ^ , ^ a 



{dQ^CTq)-L + , n Ql2. (18) 



Here the mass coefficient m2 is defined in (^. Note that without the square 
root factors in (||) the numerical coefficients in the last two terms in the first 
bracket on the r.h.s. would be — tt and 1 instead of — 37r/2 and 1/2 respectively. 

For light quarkonia the static gluon exchange interaction (|l^) has only qual- 
itative value in view of the slow convergence of the asymptotic expansion in 
v/c. Moreover the quadratic spin-orbit interaction is ill-behaved in this limit. 
A more realistic version of this interaction is obtained by replacing the main 1 /r 
(and the accompanying 5 function term) by the corresponding unapproximated 
form: 

3 r ^ 3 r ■ ^ ' 

Here the function /o(?') is defined as 

2 [°° ,^ su\(kr) M m f M + m\ ,,2, 
K eg Cg \ eg -I- / 



TT 
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The factors eg and Cq are defined as 



1,2 



(21) 



Here as(fc^) is the running coupling constant of QCD, for which we use the 
screened parametrization ||27| 



1277 



27 ln[(fc2 +4m2)/A2]' 



(22) 



The corresponding modification of the p"^ /r term in ( [l8| ) is obtained by the 
replacement 



P' 



4 

3"^2MTOr 



3 r 



(23) 



where 
/2(r) 



dfc- 



sin(fcr) /eg + A/\ f eg + 



2eQ 



265 



A/ + m 



+ Af ) (e^ + m) 2eQe^ (eg + A/)2 



Af(eQ + M) 



2p2 



+ 



1 - 



m{eq + m) 



2el 



(24) 



It is necessary to treat this term exactly in the solution of the wave equation and 
not as a first order perturbation. The unapproximated version of the spin-spin 
component of the gluon exchange interaction is introduced by the replacement 
of the delta function as 



1 



dkk sin(fcr" 



/ Af + m \ Mm 



\eQ + eqJ egCg 



(25) 



The unapproximated forms for the tensor and spin-orbit interactions in (18) are 

f M + m\ as{k'^) 



2 

Vg{T) = ^Si2 dke j2 



{kr) 



\eQ + eqJ egcq 



(26) 



Vg{LS) 



-S-L 



1 



{Cq + m)2 



dkkji{kr) 



M + m\ eQ + M eq + : 



1 - 



+ M){eq + m) 



4(eg + Af ) 
' as(k^) 



(eg + A//)2 



eg 
1 - 



4(e,- + m)2 



(27) 



The integrals in the spin-spin and tensor interactions above are strongly mod- 
erated when the nmning coupling strength as{k'^) is taken into account. The 
quadratic spin-orbit interaction component in (|l^) may be similarly regulated 
by employment of the unapproximated form 



dkk^j2{kr)- 



as{k'^ 



/ M + m 



^eqicQ + M){eq + m) \eQ 



(28) 



This form of the quadratic spin-orbit interaction has finite matrix elements. 
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The parametrization ( p2| ) takes the long distance screening of the quark- 
gluon couphng into account through the gluon mass parameter rUg. For this 
we use the value rug = 240 MeV while for the confinement scale parameter Aq 
we use the value 280 MeV. These values are chosen by a fit to the heavy-light 
meson spectra below. They are similar to those used in refs. |28[ |. 

The employment of the screened running quark-gluon coupling strength 
along with the relativistic corrections to the gluon exchange interaction leads 
to significant dampening of the gluon exchange interaction at short range. This 
is shown in Fig. |l|, where the unapproximated values for /o('') as given by the 
expression ( p^ ) are plotted for the case of D and B mesons with M — 1.58 GeV 
and 4.825 GeV respectively, and m — 450 MeV. The result reveals that the 
approximation fo{r) — as — constant is inadequate. The parametrization ( p2|) 
gives the value ~ 0.4 for as at the charmonium scale, which is consistent with 
the values extracted by lattice methods in refs. In the numerical re- 

sults for the spectra of the heavy-light mesons below, it proved essential to use 
the unapproximated expressions for the functions fo{r) and /2(?')j the reason 
being that the small mass of the light quarks render the static approximations 
misleading. 




3.3 The instanton induced interaction 



Numerical QCD lattice calculations indicate that the instanton liquid model 
provides an appropriate description of the dynamics in the low energy regime 
[H ^ . The instantons induce a pointlike interaction between heavy and light 
quarks, which contributes significantly to the hyperfine interaction in heavy-light 
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mesons ]l9| . 

The instanton induced interaction for Qq systems contains a color dependent 
component, which may be the cause of the splitting between light pseudoscalar 
and vector mesons, for which the color magnetic interaction is screened. The 
form of the instanton induced interaction, which is appropriate for heavy-light 
mesons has been derived in ref. [Q, who considered the large Nc limit, where 
the matrix element of the color exchange operator is (A^ • A^) /4 = —Nc/2. The 
spin-dependent component of this interaction couples as a T invariant 
which implies that that it may be viewed as a vector meson exchange like inter- 
action with anomalous couplings to quarks. 

If the large limit is not invoked, the instanton induced interaction for the 
Qq system takes the form [|l9| : 



Hr 



dQ-dgXQ-XgS<^^\r) (29) 



V 2nN, 
1 /AAf^P'"Am,- 



4 \ 2nNc 

Here n represents the instanton density, which typically is assigned the value 
1 fm"''. The mass parameters AAfg, Arriq and AMq™ represent the mass shifts 
of the heavy (Q) and light (q) quarks, which are caused by the instanton induced 
interaction. The numerical values for these mass shifts were obtained as AAIq= 
86 MeV, AM^"'= 3 MeV and Am,=420 MeV respectively in ref. [|9|. 

The expression ( |29| ) represents the approximate effective instanton induced 
interaction that obtains when one of the quarks is much heavier than the other 
one, which is appropriate here. The Lorcntz structure of the main term in the 
interaction is a scalar, whereas the spin dependent term has a tensor coupling 
structure. 

In the large color limit the sign of the spin-independent term in ( p^ ) is 
negative. In that limit the instanton induced interaction ( p9[ ) plays a role akin to 
that of gluon exchange, with matrix elements of similar magnitude p^ . Without 
invocation of that limit, the spin-dependent term of the interaction ( p9[ ) is not 
sufficiently strong to explain e.g. all of the difference between the D* and D 
meson masses. As below it is found that the screened gluon exchange interaction 
described above provides sufficient splitting strength, the numerical results do 
not include the interaction (|29|). If the gluon mass parameter were increased, 
it would be necessary to include the additional spin-dependence given by the 
instanton induced interaction, however. 



4 The spectra of charmonium and the heavy- 
Ught mesons 



4.1 The charmonium spectrum 

The mass of the charm quark was fixed by a fit to the D and Ds meson spectra, 
and then tested by direct computation of the cc spectrum, without further 
modification of the charm quark mass. The spectrum of charmonium is fairly 
well known experimentally. 
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The interaction model consists of the all the local terms in the scalar lin- 
ear confining interaction (^),(|8|) and the gluon exchange interaction, with full 
account of the relativistic corrections and the running coupling strength as de- 
scribed in section (^) above. 

The calculated spectrum is compared to the experimental spectrum in Ta- 
ble ^ The parameters in the interaction model are listed in Table 0. The value 
for the charm quark mass (1580 MeV) was chosen to be the same as that used 
for the D and Dg-mesons below. With the parameter values listed in Table 
|l| wc are able to accurately reproduce the experimentally determined J /'^-rjc 
splitting, but at the price of 30-50 MeV overpredictions of the excited states. 
This effect is similar to that noted in 0. Better agreement with experiment 
can be achieved by lowering the confining string constant to 960 MeV/fm and 
raising the gluon mass rrig by 20 MeV to 260 MeV. In this case the excited states 
agree fairly well with experiment, while the J/'^-rjc splitting is underpredicted 
by 15 MeV. A lower charm quark mass for the charmonium system would also 
improve the spectra, but as the charmonium spectrum is calculated primarily 
for testing purposes we keep the charm quark mass equal to that used for the 
D and Ds-mesons. We thus conclude that our test results for charmonium are 
similar to the spectrum obtained in [Q. 





cc 


Ds 


D 


c 


1120 MeV/fm 


1120 MeV/fm 


1120 McV/fm 


b 


-50 MeV 


-260 MeV 


-320 MeV 


Ao 


280 MeV 


280 MeV 


280 MeV 


rUg 


240 MeV 


240 MeV 


240 MeV 


rUc 


1580 MeV 


1580 MeV 


1580 MeV 


rUs 




560 MeV 










450 MeV 



Table 1: Model parameters used for the charmonium, D and Ds meson spectra. 



The quality of the spectrum is similar to that obtained with an effective 
interaction constructed by means of lattice methods ||l| , as well as by completely 
nonrelativistic phenomenology [ p9| [T^ . The numerical values for the calculated 
energies of the cc are listed in Table l2l along with the experimental values Jsof . 





Predicted States (MeV) 


Experimental states (MeV) 




2975 


2979.8 ± 2.1 


l^^l 


3088 


3096.88 ± 0.04 


2^50 


3682 




2^Si 


3736 


3686.00 ± 0.09 


I'Pi 


3518 




l^Po 


3450 


3417.3 ± 2.8 


l3Pl 


3519 


3510.53 ± 0.12 


13P2 


3580 


3556.17 ± 0.13 



Table 2: Calculated and experimental charmonium states. 
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4.2 The D and Dg -meson spectra 



In the calculation of the spectra of the D and Dg mesons, the constituent masses 
of the light and strange quarks are treated as phenomenological parameters to be 
fitted against the known splittings in the spectrum. While a reasonable overall 
description of the D meson spectrum is readily obtainable with the interaction 
models described above, a quantitative description of e.g. the spin-orbit splitting 
of the P-shell Dj {j = 0, 1, 2) D meson resonances demands a real effort. 

The main difhculty with the spin orbit splittings may be ascribed to the 
small mass of the light quarks, which makes the matrix elements of the spin- 
orbit components of both the confining and hyperfine interactions large. That is 
the reason for taking into account the relativistic damping of these interaction 
components that is described in section ^. 



13^1 

2' So 
2^ Si 
3' So 
3^ Si 
'i'So 
4^Si 



l^Pi 

2^Pi 

2'Po 
2^Pi 

23P2 

3iPi 
S^Pi 



l^Di 

2^D2 
2^Di 
2'^D2 
2^D3 



Mu d=450 MeV 



1874 
2006 
2540 
2601 
2904 
2947 
3175 
3208 



2389 
2341 
2407 
2477 
2792 
2758 
2802 
2860 
3082 
3050 
3085 
3142 



2689 
2750 
2727 
2688 
2997 
3052 
3029 
2999 



Exp(j:)") Exp(j:)^) 



1864.6±0.5 
2006.7±0.5 

2637 ? 



1869.3±0.5 
2010.0±0.5 

2637 ? 



2422. 2±1. 8 
2458.9±2.0 



2459 ± 4 



Table 3: Calculated and experimental D-meson states. 
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The D Spectrum (State energies in MeV) 



Experimental States + 
IVI = 1 580 MeV, m = 450 MeV x 



'D, ■'D, ■'D, 



Figure 2: Calculated and experimental D-meson states. 



The spectrum of the D meson that is obtained with the scalar linear confining 
interaction and with the quark masses rriu = nid = 450 MeV and rric — 1580 
MeV is shown in Fig. ^ The light quark mass here is about 100 MeV larger 
than the typical values employed in nonrelativistic phenomenology. Reducing 
this mass further in the calculation would lead to an unrealistically large spin- 
orbit splitting of the P— states, while giving an unrealistically small IS — > 2S 
spacing. The interaction parameters used in the calculation are listed in Table |l]. 

The calculated energies of the Z?— meson states are listed in Table ^ along 
with the known empirical energies |^o| . 

The calculated D meson spectrum has the correct D* — D ground state 
splitting and only slightly underpredicts the possible D meson excitation 
at 2637 MeV. The P-shell states around 2400 MeV are also well reproduced. 

The calculated energies of the Ds meson states are shown in Fig. |^ and are 
listed in Table ^. As in the case of the D meson a satisfactory description of the 
still very incompletely known experimental spectrum is only achievable with a 
fairly large value for the constituent mass of the s quark (tos = 560 MeV). 

The overall structure of the Dg meson spectrum in Fig. ^ is similar to that of 
the D meson spectrum in Fig. ^j. The ground state splitting is given correctly, 
as are the known P-state excitations. 

The calculated D meson spectrum is similar to that obtained by the Gross 
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reduction of the Bethe-Salpeter equation in ref.|^. The shell spacing at in- 
creasing energy of the latter is somewhat wider, which may reflect a difference 
between the different quasipotential reductions, but also the employment of 
static potentials in ref.P]. 





Ms=560 MeV 


Experimental states (MeV) 




1975 


1968 5 ± 6 




2108 


2112.4 ± 0.7 




2659 




2^ Si 


2722 




3^50 


3044 




3351 


3087 




4:^ So 


3331 




4^Si 


3364 




l^Pi 

i 


2502 




l^Pn 
^ u 


2455 






2522 


2535.35 ± 0.34 ± 0.5 




2586 


2573.5 ± 1.7 


2iPi 


2928 






2901 




23Pi 


2942 




23P2 


2988 




31P1 


3234 




33^0 


3214 




33Pi 


3244 




33P2 


3283 






2838 






2845 




1=^P'2 


2856 






2857 




2^D2 


3144 




2^Di 


3172 




2^D2 


3167 




2^D3 


3157 





Table 4: Calculated and experimental Ds-meson states. 
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The Dg Spectrum (State energies in MeV) 



Experimental States + 
M = 1 580 MeV, m = 560 MeV x 



X X 



'So ^Si Vi ¥(, ^P, ^Pj ^D, ^D; 



Figure 3: Calculated and experimental Dg meson states. 



4.3 The B and Bg -meson spectra 

The empirical knowledge of the spectra of the B and Bg mesons is still very 
incomplete. Only the ground state pseudoscalar and vector meson state energies 
are known with certainty. In addition one orbital excitation of the B meson is 
known at 5700 MeV, which presumably belongs to the P— shell. 

The only additional parameter required to calculate the spectra of the bot- 
tom mesons with the interaction models described above is the b quark mass. 
For this we choose the value mb=4825 MeV, with which value the pseudoscalar- 
vector splitting of the bottom and strange bottom mesons is obtained correctly. 

The calculated B and Bs meson spectra are shown in Figs. ^ and ||. The 
energies of these states are also listed in Tables |^ and|7[ along with the empirical 
values that are taken from ref.[0. The orbital excitation of the B meson at 
5700 MeV is obtained correctly, under the assumption that it corresponds to a 
j = I P-shell state. 

The overall features of the calculated bottom meson spectra are similar to 
those obtained with the Gross reduction [Q of the Bethe-Salpeter equation in 
ref.Pl, although as in the case of the charm mesons, we obtain somewhat smaller 
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shell spacings at higher excitation. 





Bs 


B 


c 


1120 MeV/fm 


1120 MeV/fm 


b 


-185 MeV 


-250 MeV 


Ao 


280 MeV 


280 MeV 


rUg 


240 MeV 


240 McV 


mb 


4825 MeV 


4825 MeV 


rus 


560 MeV 




mu,d 




450 MeV 



Table 5: Parameter values used in the calculation of the Bg and S-meson 
spectra. 



State 



M„ d=450 MeV 



Exp(E") Exp(g±) 



13^1 

23^1 

31^0 
43^1 



l3Po 

2^ Pi 
23Po 

23Pi 

2^P2 

33Po 

33Pi 
S^Po 



5277 
5325 
5822 
5848 
6117 
6136 
6335 
6351 



5279. 2±1.^ 
5324.9±1.^ 



5278.9±1.8 



13D2 

2^D2 
2^Di 
2^D2 

23£>3 



5686 
5678 
5699 
5704 
6022 
6010 
6028 
6040 
6259 
6242 
6260 
6277 



5697±9 



5920 
6005 
5955 
5871 
6179 
6248 
6207 
6140 



Table 6: Calculated and experimental B-meson states 
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Ms=560 MeV 


Experimental Bs states (MeV) 


1^5*0 


5366 


5369.3 ± 2.0 




5417 


5416.3 ± 3.3 




5939 




2^ Si 


5966 




3^ So 


6254 




3^ Si 


6274 




4:'^ So 


6487 




4^ Si 


6504 




l^Pi 


5795 






5781 






5805 






5815 




2^ Pi 


6153 




2^Pq 


6143 




2^ Pi 


6160 




2^P2 


6170 




3^ Pi 


6406 




3^Po 


6396 




3^ Pi 


6411 






6421 






6043 




l^Di 


6094 






6067 






6016 




2^D2 


6320 




2^Di 


6362 




2^D2 


6339 






6298 





Table 7: Calculated and experimental Sg-meson states. 
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The B Spectrum (State energies in MeV) 











1 






Experimentai States + 
= 4825 MeV, m = 450 iWeV x 


6400 


X 


X 


















X 


X 


X 


X 


X 


6200 


X 


X 










X 

X 

X 


6000 






X 


X 


X 


X 


X 

X 

X 


5800 


X 


X 


X 


X 


X 


X 


X 


5600 
















5400 




X 












5200 








1 








\ 




'P1 


3p 




3p 


^2 ^D, ^Dg ^Dj 



Figure 4: Calculated and experimental B-meson spectra. 
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The Bs Spectrum (State energies in MeV) 



Experimentai States + 
M = 4825 MeV, m = 560 iWeV x 



^Sq ^S, 'p, ^Pq ^Pi ^P2 ^Dg ^D, ^Dg ^Dj 



Figure 5: Calculated and experimental -Bg-meson spectra. 
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5 Quark current operators 



5.1 Single quark currents 

In order to calculate the Ml transitions of the Qq systems we employ the full 
Dirac structure of the single quark current operators. This was shown to be 
necessary in ref.|l^ for a satisfactory description of the corresponding Ml de- 
cays of heavy quarkonia. The current density operators of the constituent light 
quarks and the c and b quarks are then 

{p'\m\p) = tQ[uip')ju{p)]. (30) 

Here Q is the quark charge operator, which takes the values ±2e/3 and ±e/3. 
With canonical boosts for the quark spinors u, u, this current operator leads to 
the magnetic moment operator 

where v is the quark velocity: v — {p' + p)/2mQ, rup is the proton mass and 
is the nuclear magneton. 

For the mesons the spin flip part of the combination of the magnetic 
moments of the charm quark and the light antiquarks is 



1 + 



^J■N■ 



(32) 



For the Z?°, D'^ mesons the corresponding operator is 




mu y/1 + vl 



Here m is the mass of the light antiquark, and Vc and Vq denote the velocity 
operators of the charm quark and the light antiquark respectively. The spin- 
flip magnetic moment operator for the Df^ meson may be obtained from the 
corresponding expression for the mesons by replacing ctj, and in 



52|) by CTj, rUs and Vg respectively. 

The spin-flip magnetic moment operators for the bottom mesons may be 
constructed by similar replacements. For the mesons these operators are 
obtained from those for the ( |3^ ) by replacement of the c quark mass and 
velocity by the corresponding u quark mass and velocities and by replacement 
of the d quark mass and velocity by the b quark mass and velocity. 

The spin-flip magnetic moment operator for the and mesons is ob- 
tained from that for the neutral D mesons (B3[) by multiplying the operator by 



20 



— 1/2 and by replacing the c and u quark masses and velocities by the corre- 
sponding masses and velocities for the d and b quark. A further replacement of 
the d quark mass and velocity by the corresponding s quark mass and velocity 
yields the spin-flip magnetic moment operator for the Bg mesons. 

Note that when matrix elements of these magnetic moment operators are 
evaluated with wave functions that are solutions to the Blankenbecler-Sugar 
equation (|^), these operators should be multiplied by the factor 

fBS = ^= (34) 

MJl + v'^+mJl + vl 

5.2 Two-quark currents 

Both the scalar confining interaction and the hyperfine interaction may excite 
virtual quark-antiquark states, which are deexcited by the electromagnetic field. 
This process is illustrated by the Feynman diagrams in Fig.^ . These operators 
may be derived by the methods described in ref. p6| . 




Figure 6: Exchange current operators associated with the effective scalar con- 
fining and gluon exchange interactions, with intermediate virtual qq excitations. 

For the system the two-quark current that is implied by the scalar con- 
fining interaction to lowest order in v/c takes the form 

UC) = T*ecr^ (A _ J_\ (^^ „ g^) x q. (35) 

12 \ TOJ TOj / 

This current operator may be viewed as a direct renormalization of the corre- 
sponding spin-flip combination of the single c and d current operators: 

m ( 2 1 \ 

J,, = ±^e^ <fa) X q. (36) 
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The two-quark magnetic moment operator that is imphed by this current 
operator is then accordingly 

m f 2 1 1 

M2(C) Tcr— ^ {— 2 r (f^c - crj) fiN- (37) 

For the D^,D^ mesons the corresponding two-quark magnetic moment is 
then 

M2(C) = Tcr^(^-^|(^c-afi)MAr. (38) 

The exchange magnetic moment operator for the Df- mesons is obtained from 
that for the mesons (^7|) , by replacing ctj and m j by CTj and respectively. 

The general rule is that the exchange current that is associated with the 
confining interaction has opposite sign to the corresponding combination of 
single quark current operators, and may be obtained directly from the former 
by multiplication by the factor — cr and by squaring the quark masses in the 
denominators. The corresponding operators for the bottom mesons are readily 
constructed with the help of this rule. Note that the constant 6, which is 
subtracted from the Hamiltonian, is viewed as an approximation to the velocity 
dependent term — 3/2p^/m2 in the confining interaction (^. As that term may 
be viewed as the origin for the two-quark interaction current that is associated 
with the scalar confining interaction through the continuity equation, it should 
not be subtracted from the cr term in the two-body current. 

The gluon exchange interaction (|l^) also implies a two-quark exchange cur- 
rent operator. The matrix elements of this two-body current are however sig- 
nificantly weaker than those of the exchange current operator that is associated 
with the scalar confining operator jl^ . Because of the cancellations between the 
matrix elements of the single quark current operators and those of the exchange 
current operator that is associated with the confining interaction, the gluon ex- 
change current operator does give a significant contribution to the calculated 
Ml transition rates between the ground state heavy-light mesons. 

For the mesons this two-quark current, to lowest order in the inverse 
quark masses, has the form 



UG) \D' 



'Pd+Pd 




(39) 



Here the momentum operators /cj and denote the fractional momenta 
imparted to the d and c (or d and c) quarks respectively. 

The spin-fiip component of the magnetic moment operator, that is obtained 
from this current operator, is then 



UG) [D^] 



^ .fojr) nip I 
27 r y 



(40) 
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In the case of the and mesons the gluon exchange current operator 
takes the form 



±e- 



IGtt f as{kl) 



27 



Pt, 



Pt, 



CTu 



mcniu 



mcrriu 
X ku 



Pc 



rricTnu 



+ i 



The corresponding spin-flip magnetic moment operator is 



2/0 (r) mp 
27 r 



1 ^ \(- - \ 

K T \Cc — CTfij MW- 



.(41) 



(42) 



The gluon exchange spin-flip magnetic moment for the mesons is ob- 
tained from that for the mesons ( ^0| ) by replacing ctj and toj by CTj and mj 
respectively. 

The gluon exchange spin-flip magnetic moment operator for the -B" and 
mesons is obtained from the corresponding expressions ( |4^ ) for the 13° and 

mesons by replacing the c quark mass and spin operator by the h quark 
mass and spin operator. Similarly by replacing the c quark mass and spin 
operator in the expressions (^) for the D"*" and D~ gluon exchange spin-flip 
magnetic moment operator by the h quark mass and spin operator one obtains 
the corresponding operator for the and JS" mesons. Finally the expressions 
for the gluon exchange spin-flip magnetic moment for the -B^ and -B^ mesons 
are obtained from those for the and D~ mesons by the replacement of the 
c quark mass and spin operator by the h quark mass and spin operator. 

The exchange magnetic moments that are associated with the conflning and 
gluon exchange interactions above represent relativistic corrections of second 
order in vjc^ as is implied by their origin as pair excitation currents (Fig. 
They thus appear in the same order as the relativistic corrections to the single 
quark magnetic moment operators. As relativistic corrections to the exchange 
current operators are of quartic order in vj c we have not considered them here, 
in order to be consistent with the treatment of the conflning interaction without 
quartic terms in w/c, besides the numerically insigniflcant quadratic spin-orbit 
interaction. 



6 Ml transitions 



6.1 Ml transition matrix elements for D mesons 

The spin-flip component of all the magnetic moment operators derived above 
may be written in the general form 

ft^eM ((TQ - aq) . (43) 

Here Ai is the matrix element of the orbital part of the magnetic moment 
operator between the initial and final meson states. The transition width for 
Ml transitions of the type Qq{J = 1) ^ Q<i{J — 0)7 may be expressed as 

r = yae™^A^^g^ (44) 
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Here aem is the fine structure constant and q is the photon momentum in the 
laboratory frame. The masses of the initial and final meson states are denoted 
Mi and Mf respectively. 

The matrix element Ai is defined as the coefficient of the spin-flip operator 
Bc—dq in the expression for the magnetic moment operator. It is formed as a sum 
of matrix elements of single quark and exchange magnetic moment operators. 
In the case of the mesons, the matrix element between 5-states of the single 
quark operator (|3^ ) may be written in the form 



Mia{D^) 




drr'^ipf{r') 



(45) 



jo (r'y/p^ +pqz/2 + q'^/16j jo (^ry/p'^ - pqz/2 + q'^/ldj ipi{r) 

Here Vc — p/rric^ = p/m^, and the factor fssip) is the "minimal relativ- 
ity" factor (34). The expression for corresponding matrix elements in the case 
of the 1?° and systems is readily inferred from eqn ( p3| ) by comparison. Here 
ipi{r) and ff{r') denote the initial and final orbital wavefunctions respectively. 

It is instructive to compare these expressions to those of the static quark 
model (uc, Vj; 0). In the case of the mesons that expression is simply 



(46) 



whereas for the D , systems it is 



Mia{D'',D') = ±i r d\^)\r) (— + — ) JO (%) V^{r). (47) 

D Jo xTTlc nifi / \ A ' 

The expressions (^ ) and ( ^6|) apply to mesons with the substitutions 
nig and Vg. 

The matrix elements of the exchange magnetic moment that is associated 
with the confining interaction (||) is in the case of the mesons 

M.iC) = ^Icfd^rr^^r) (j, - i,) ,o (f ) ..(0- (48) 

This expression applies to the mesons as well, with the substitution 
mg. In the case of the Ug, Dq mesons the corresponding expression is 

M,(C) . ^cf d=.,v;(r) [± + Jj) i. (f ) ,.(r), (49) 

The matrix element of the gluon exchange spin-flip magnetic moment oper- 
ator is in the case of the mesons (cf.(^)): 
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In the case of the D'-' , mesons this matrix element is 



The rule for constructing the corresponding matrix element expressions for 
the other heavy-light meson charge states may be inferred from these expres- 
sions. 



6.2 Ml transitions of charm mesons 



At present the only experimental information on the Ml transition rates of D- 



meson states are the Ml branching ratios for £'*(2010)^ 
for i:>*(2007)'' 



■7: 



1.1 



+2.1 



% and 



D^r- 38.1 ±2.9% [30 



As the upper limit for the total width 
of the former transition is 0.131 MeV and for the latter 2.1 MeV, the upper 
limits for the Ml decay widths are accordingly 1.44 keV for Z?*(2010)^ D^j 
and 800 keV for D*(2007)° -> D°-f. These upper hmits suggest that the Ml 
decay widths should be of the order 1 keV. 





NRIA 


RIA 


RIA + cf. 


RIA + cf. + ex. 






20.9 keV 


6.52 keV 


2.21 keV 


1.25 keV 


D*'° - 




805 eV 


51.1 keV 


39.3 keV 


44.5 keV 


D*'° - 




1.88 keV 


631 eV 


2.48 keV 


2.24 keV 






3.02 keV 


35.7 keV 


31.6 keV 


34.2 keV 



Table 8: The Ml decay widths of the Dq (cu) mesons 





NRIA 


RIA 


RIA + cf. 


RIA + cf. -1- ex. 






574 eV 


10.8 eV 


2.09 keV 


1.10 keV 


£)*'± _ 




23.9 eV 


2.77 keV 


12.5 keV 


15.7 keV 


£)*'± ^ 




52.6 eV 


2.84 eV 


748 eV 


607 eV 


z?'±- 




85.3 eV 


1.63 keV 


10.2 keV 


11.9 keV 



Table 9: The Ml decay widths of the D± [cd) mesons 



This is indeed what we find by explicit calculation here. The calculated 
widths for the transitions D*(2010)± D^-f and D*(2007)° L»°7, as well 
as for the transitions between their respective excited states are listed in Ta- 
bles ^ and H respectively. The first column gives the nonrelativistic impulse 
approximation results (NRIA), which are obtained with static quark current 
operators. The experience from the corresponding decays in charmonium sug- 
gests that these are not realistic. The results obtained in the relativistic impulse 
approximation, with the full Dirac magnetic moment operators are listed in the 
second column (RIA). These values should also, on the basis of the experience 
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NRIA 


RIA 


RIA + cf. 


RIA + cf. + ex. 






178 eV 


7.10 mcV 


736 eV 


337 eV 


D* - 




7.88 eV 


1.22 keV 


3.26 keV 


4.47 keV 


D* - 




17.1 eV 


246 meV 


258 eV 


200 eV 


D's- 




22.4 eV 


654 eV 


2.53 keV 


3.13 keV 



Table 10: The Ml decay widths of the Dg (cs) mesons 



with heavy quarkonia, be unreahstic. FinaUy the column RIA + cf. give the 
results that are obtained when the exchange current contribution that is as- 
sociated with the scalar confining interaction is taken into account. The last 
column (RIA + cf. + ex.) also takes into account the contribution of the gluon 
exchange current. The contribution of the latter is most significant for the Ml 
transitions between the ground states. 

In Tables ^ and |^ the Ml transition between the excited D* and D mesons 
and the ground states are given as well. In the case of the forbidden transitions 
between the orbitally excites states the two-body current contributions are dom- 
inant. The large difference between the calculated Ml transition rates of the 
charged and neutral D-mesons is due to the very different quark mass depen- 
dence of the associated spin-flip magnetic moment operators. In the derivation 
of these results we determined the value of the photon momentum q from the 
differences between the empirical values for the D meson states as given in 
ref.[0. The difference is only notable in the case of the excited D' and D * 
mesons, the latter of which is underpredicted by about 35 MeV. 

We have also calculated the Ml transition rates for J/-0 — > 77c7 and 'ip' rjcj 
with the present framework and parameters. In this case the gluon exchange 
current does not contribute at all because of the equal masses of the charm 
quarks and antiquarks. The empirical width for the J/ip — > ric"f decay is known 
to be 1.14 ± 0.39 keV. For this we find the value 0.75 keV when the two-body 
contribution is included, while the relativistic impulse approximation value is 
1.47 keV. The static quark model prediction is 2.51 keV, which amounts to an 
overprediction of about a factor 2. The empirical width of the forbidden Ml 
transition — > r/c7 is 0.78 ± 0.24 keV. For this we obtain the value 1.49 keV 
when the two-body current is taken into account. This is much closer to the 
empirical value than the value 0.17 keV given by the static quark model or the 
value 8.48 keV that is obtained in the relativistic impulse approximation. 

The Ml calculated transitions between the D* and Dg mesons are listed in 
Table |lO|. These transition rates should a priori be expected to be qualitatively 
similar to those of the D*^ and mesons. This is borne out by the calcu- 
lation, although the numerical values show marked differences. This is another 
illustration of the sensitivity of the calculated Ml transition rates to the pa- 
rameters of the model and in this case to the difference between the constituent 
masses of the s and d quarks. 

6.3 Ml transitions of bottom mesons 

The calculated Ml transition widths of the bottom mesons are listed in Tables ^ 
and ^ The experimental decay widths for these transitions are yet to be 
measured. The large masses of the bottom mesons make these transition widths 
much smaller than those of the corresponding transitions of charm mesons. 
The overall features of these results are however similar to those of the 
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charm mesons. The nonrelativistic impulse approximation leads to considerable 
overestimates. In the case of the forbidden transitions the two-body exchange 
currents are again dominant. 

The calculated Ml transition rates for the strange bottom mesons are listed 
in Table 





NRIA 


RIA 


RIA + cf. 


RIA + cf. + ex. 


B*" - 


^B^ 


152 eV 


51.9 eV 


25.1 cV 


9.55 eV 


B*'" - 




149 meV 


5.39 keV 


9.73 keV 


12.2 kcV 


B*'° - 




27.7 eV 


9.59 eV 


48.3 eV 


40.6 eV 


fiO' - 


, 5*0 


34.6 eV 


4.31 keV 


7.95 keV 


9.72 keV 



Table 11: The Ml decay widths of the Bq (bd) mesons 





NRIA 


RIA 


RIA + cf. 


RIA + cf. + ex. 


B*± - 




462 eV 


132 eV 


160 eV 


67.4 eV 


B*'± - 




455 meV 


20.1 keV 


39.2 keV 


50.9 keV 


B*'± - 




84.3 eV 


24.6 eV 


221 eV 


183 eV 


b'± - 




105 eV 


15.8 keV 


32.4 keV 


40.7 keV 



Table 12: The Ml decay widths of the B^ {bu) mesons 





NRIA 


RIA 


RIA + cf. 


RIA + cf. + ex. 


B*- 


^Bs 


116 eV 


46.4 eV 


945 meV 


148 meV 


B*' - 


Bs 


2.39 eV 


3.87 keV 


2.29 keV 


3.22 keV 


B*' - 




20.9 eV 


8.26 eV 


12.4 eV 


9.76 eV 




.b: 


10.6 eV 


3.12 keV 


1.84 keV 


2.50 keV 



Table 13: The Ml decay widths of the Bg {bs) mesons 



7 Discussion 



The main result of the present study is that it is possible to achieve a satisfac- 
tory description of the heavy-light meson spectra, as well as credible predictions 
for their radiative spin-flip transitions with the framework of the Blankenbecler- 
Sugar equation with a scalar confining interaction. This conclusion differs from 
that in ref. |Tl| . The calculated results are however very sensitive to the param- 
eters in the model, and to the relativistic damping of the hyperfine interaction. 
This suggests that that model for the hyperfine interaction (relativistic gluon 
exchange) may not provide a sufficiently complete description of the hyperfine 
interaction. It may be conjectured that screening the gluon exchange interaction 
at shorter distances and compensating with the instanton induced interaction 
may lead to less parameter sensitive results. 
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Both the Blankenbecler-Sugar and Gross quasipotential reductions of the 
Bethe-Salpeter equation lead to two-body exchange current operators, which 
connect to negative energy intermediate states. If the confining interaction cou- 
ples as a scalar to the quarks, the accompanying two-quark magnetic moment 
is a pure spin-flip operator for mesons with equal quark and antiquark masses. 
The contribution of this operator is essential for obtaining realistic Ml transition 
rates for heavy quarkonia |l6| . In the equal mass case the analog two-quark op- 
erator that arises with a vector confining interaction has no spin-flip component 

For mesons with unequal quark and antiquark masses the gluon exchange 
current operator has a spin- flip component. It therefore also contributes to the 
Ml decay rates of the heavy-light mesons. Although this contribution is much 
smaller than that of the confining interaction, it is significant in the case of the 
Ml transitions between the ground state vector and pseudoscalar heavy-light 
mesons, because of the partial cancellation between the matrix elements of the 
single quark current operators and those of the exchange current associated with 
the confining interaction. 

Strong conclusions concerning the Ml transition rates here have to await 
experimental determination of the total width of the heavy-flavor vector mesons. 
As the relative branching ratios for Ml decays of these mesons are known, the 
key missing empirical information is their total widths for strong decay. 
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